Mobility edge phenomenon in a Hubbard chain: A mean field study 
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We show that a tight-binding one-dimensional chain composed of interacting and non-interacting 
atomic sites can exhibit multiple mobility edges at different values of carrier energy in presence of 
external electric field. Within a mean field Hartree-Fock approximation we numerically calculate 
two-terminal transport by using Green's function formalism. Several cases are analyzed depending 
on the arrangements of interacting and non-interacting atoms in the chain. The analysis may be 
helpful in designing mesoscale switching devices. 
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I. INTRODUCTION 

Electronic localization phenomena in onc-dimcnsional 
(ID) quantum systems have long been a central prob- 
lem in condensed matter physics. It is well established 
that in infinite ID systems with random site potentials, 
irrespective of the strength of randomness, all the energy 
cigenstates are exponentially localized^. Apart from this 
Anderson type localization another kind of localization 
known as Wannier-Stark localization is also observed in 
ID materials, even in absence of any disorder, when the 
system is subjected to an external electric field 2 -. For 
both cases i.e., infinite ID systems with random site 
potentials and ID chains in presence of external elec- 
tric field, one never encounters mobility edges separating 
the localized energy eigenstates from the extended ones, 
since all eigenstates are localized. However, there are 
some classes of ID systems such as, correlated disordered 
models, quasi-periodic Aubry- Andre model where several 
classic features of mobility edges at some specific values of 
energy are obtained^— . Although the existence of such 
mobility edges in one- or two-dimensional systems has 
been described by several groups^—, a comprehensive 
study of this phenomenon is still lacking, particularly in 
the presence of electron-electron interaction. Still open 
fundamental questions are whether some special features 
exist in disordered ID systems, or in the response of ID 
systems to an externally applied electric field even in the 
presence of electron-electron interaction. 

In the present article we investigate two-terminal elec- 
tron transport through a ID mesoscopic chain composed 
of interacting and non-interacting atomic sites in pres- 
ence of external electric field. Although some works have 
been done in such superlattice structures^—, the anal- 
ogous representation of metallic multilayered structures 
which exhibit several novel features^— , no rigorous ef- 
fort has been made so far, to the best of our knowledge, 
to unravel the effect of the interplay of electron-electron 
interaction and an imposed external electric field on elec- 
tron transport in such systems. Here we show that a 
traditional ID lattice with electron-electron interaction, 
evaluated at the Hartrcc-Fock (HF) mean field (MF) 



level, is characterized by a mobility edge behavior at fi- 
nite bias voltage. Furthermore, a superlattice structure 
comprising sites on which electron-electron interactions 
are expressed differently (some sites are interacting and 
some sites are non-interacting) is characterized by multi- 
ple occurrence of mobility edges at several values of the 
carrier energy. The applicability of mean field approxi- 
mation in such superlattice geometries has already been 
reported in a recent work 2 ^. 



II. MODEL AND CALCULATION 

We adopt a tight-binding (TB) framework to describe 
the model quantum system and numerically calculate 
two-terminal transport within a mean field Hartree-Fock 
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FIG. 1: (Color online). A ID mesoscopic chain, composed 
of interacting (filled black circle) and non-interacting (filled 
green circle) atomic sites, is attached to two semi-infinite ID 
metallic electrodes, representing source and drain. 

approximation using a Green's function formalism. Sev- 
eral cases characterized by different arrangements of in- 
teracting and non-interacting atomic sites in the chain, 
are analyzed. For these models we calculate the average 
density of states (ADOS) and the two-terminal transmis- 
sion probability, and find that sharp crossovers from com- 
pletely opaque to fully or partly transmitting zones take 
place at one or more specific electron energies. This ob- 
servation suggests the possibility of controlling the trans- 
mission characteristics by gating the transmission zone, 
and using such superlattice structures as switching de- 
vices. 

Let us refer to Fig. [TJ where a ID mesoscopic chain, 
composed of non-interacting and interacting atomic sites, 
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is attached to two semi-infinite ID non-interacting source 
and drain electrodes. In the arrangement of the two dif- 
ferent atomic sites shown in Fig. Q] M (M > 1) non- 
interacting sites are placed between two interacting sites. 
Here and in what follows we make a restriction that in- 
teracting atoms are not placed successively. In a Wannier 
basis, the TB Hamiltonian for a iV-site chain reads, 

H C = Y e ™ C L C ™ + Y * [4a C 3<y + C )a C i<y 
i,a (ij),a 

+ Y UiC^ctfC^cn (1) 

i 

where, c icr (ci a ) is the creation (annihilation) operator 
of an electron at the site i with spin a (=t)4-)j t 1S the 
nearest-neighbor hopping element, ej CT is the on-site en- 
ergy of an electron at the site i of spin a and Ui is the 
strength of on-site Coulomb interaction where C7, = 
for the non-interacting sites. In presence of bias voltage 
V between the two electrodes an electric field is devel- 
oped and the site energies become voltage dependent, 
€i a = e? + £i(V), where e° is a voltage independent term. 
For the ordered chain e° is a constant independent of i 
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FIG. 2: (Color online). Variation of voltage dependent site 
energies in a ID chain with 300 lattice sites for three different 
electrostatic potential profiles when the bias voltage V is fixed 
at 0.2. 

that can be chosen zero without loss of generality, while 
for the disordered case we select it randomly from a uni- 
form "Box" distribution function in the range — W/2 to 
W/2. 

The voltage dependence of Cj(V) reflects the bare elec- 
tric field in the bias junction as well as screening due to 
longer range e-e interaction not explicitly accounted for 
in Eq. [T] In the absence of such screening the electric field 
is uniform along the chain and e»(V) = V/2 — iV/(N+l). 
Below we consider this as well as screened electric field 
profiles, examples of which are shown in Fig. [2J We will 
see that the appearance of multiple mobility edges in su- 
perlattice geometries strongly depends on the existence 
of finite bias and on the profile of the bias drop along the 
chain. 

The Hamiltonian for the non-interacting (Ui = 0) elec- 
trodes can be expressed as, 

^lead = Y e °4 c P + Y to ( c l c i + c I c p) ( 2 ) 
p <p<?> 



with site energy and nearest-neighbor intersite coupling 
eo and to, respectively. These electrodes are directly cou- 
pled to the ID chain through the lattice sites 1 and N. 
The hopping integrals between the source and chain and 
between the chain and drain are denoted by ts and td, 
respectively. 

In the generalized HF approach^—, the full Hamilto- 
nian is decoupled into its up-spin and down-spin compo- 
nents by replacing the interaction terms by their mean 
field (MF) counterparts. This redefines the on-site ener- 
gies as e' it = e, t + [/(n^) and e' 4 = e^ + U (n^) where, 

nta = cl a Ci a is the number operator. With these site en- 
ergies, the full Hamiltonian (Eq. [1]) can be written in the 
MF approximation in the decoupled form 

H MF = Y e ^t + Y f [ c lt c rt + 4t Ci t 

+ Y e W + Y 1 [ c k c u + 

' (if) 
i 

= Hc,t + Hc^-Y U ^)( n n) ( 3 ) 

i 

where, i?c7,t anc ^ correspond to the effective TB 




FIG. 3: (Color online). Transmission probability T (red color) 
and ADOS (green color) as a function of energy E for a ID 
non-interacting (Ui = V i) ordered (W = 0) chain with N = 
300 sites. The electrostatic potential profile varies linearly 
(red curve in Fig. [2}, with the total potential drop across the 
chain to be (a) V = and (b) V = 0.2. 

Hamiltonians for the up and down spin electrons, respec- 
tively The last term provides a shift in the total energy 
that depends on the mean populations of the up and 
down spin states. 

With these decoupled Hamiltonians (i?c,t an d Hc,i) 
of up and down spin electrons, we start our self consis- 
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FIG. 4: (Color online). Transmission probability T (red color) and ADOS (green color) as a function of energy E for an 
ordered (W = 0) ID chain. The left column corresponds to the case where all sites are interacting (Ui = 2), while the right 
column represents the results for a ID superlattice geometry where four non-interacting (Ui = 0) atoms are placed between 
two interacting (Ui = 2) atoms. The 1st, 2nd and 3rd rows correspond to V = 0, 0.1 and 0.2, respectively. All these results 
are shown for a linear bias drop along the chain. 



tent procedure considering initial guess values of (n^) 
and (nil). For these initial set of values of (n^) and 
(nil}, we numerically diagonalizc the up and down spin 
Hamiltonians. Then we calculate a new set of values of 
(njf) and (nq.). These steps are repeated until a self 
consistent solution is achieved. 

The converged mean field Hamiltonian is a sum of sin- 
gle electron up and down spin Hamiltonians. The trans- 
mission function is therefore a sum T(E) = '^2 cr T cr (E) 
where2£ T a = Tr [r s G r C(J T D G a c a ] . Here, G r C c and 
G c q are the retarded and advanced Green's functions, 
respectively, of the chain including the effects of the elec- 
trodes. Gc,a = (E — Hc\a — — £_d)~\ where Eg and 
YiQ are the self-energies due to coupling of the chain to 
the source and drain, respectively, while Y$ and Tr> arc 
their imaginary parts. 



III. RESULTS AND DISCUSSION 

In what follows we limit ourselves to absolute zero tem- 
perature and use the units where c = h = e = 1. For the 
numerical calculations we choose t = 1, eo = 0, to = 3 
and Ts = To = 1. The energy scale is measured in unit 
oft. 

Before addressing the central problem i.e., the possi- 



bility of getting multiple mobility edges in ID superlat- 
tice geometries, first we explore the effect of finite bias 
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FIG. 5: (Color online). Transmission probability T (red color) 
and ADOS (green color) as a function of energy E for a ID 
chain (N — 300) in absence of disorder (W = 0) with on- 
site interaction Ui = 2 and bias voltage V — 0.2 that varies 
linearly along the chain. Here we set M = 5 i.e, five non- 
interacting atoms are placed between two interacting atoms. 

on electron transport in two simple systems, one for a 
standard non-interacting chain and the other for a con- 
ventional Hubbard chain where all sites are interacting. 

In Fig. [3] wc show the variation of total transmis- 
sion probability (T) together with the average density of 
states as a function of energy E for an ordered (e° = for 
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all atomic sites i in the chain) non-interacting chain for 
two different magnitudes of the voltage bias, assuming a 
linear bias drop (uniform electric field) across the chain. 
In the absence of electric field electron transmission takes 
place throughout the energy band as clearly seen from 
the spectrum Fig. (3][a) , since in this case all the energy 
cigcnstates are extended. On the other hand, when a 
finite bias drop takes place along the chain, several en- 
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FIG. 6: (Color online). Same as Fig. [5] with M = 6. 

ergy eigenstatcs appear in the energy regions around the 
band edges for which the transmission probability is ex- 
actly zero (Fig. [3Jb)). Therefore, the chain appears in- 
sulating when Fermi energy is within the zone of zero 
transmission, while finite transmission, T ^ 0, is seen 
more towards the band centre. The sharp transition be- 
tween these regimes illustrates the existence of a mobility 
edge phenomenon under finite bias condition. For a fi- 




FIG. 7: (Color online). Transmission probability T (red color) 
and ADOS (green color) as a function of energy E for a ID 
chain (iV = 300) in presence of disorder (W — 0.5) for the 
same parameter values used in Fig. 2Jf): Ui = 2, V — 0.2 
(linear potential profile) and M = 4. 

nite bias, the localization of energy levels always starts 
from the band edges and the width of the localized energy 
zones can be controlled by the imposed electric field. Ob- 
viously, for strong enough electric field almost all the en- 
ergy levels are localized and the extended energy regions 
disappear, so that in this particular case metal-insulator 
(MI) transition will no longer be observed. This local- 
ization phenomenon in presence of an external electric 
field has already been established in the literature, but 
the central issue of our present investigation - the in- 
terplay between the Hubbard interaction strength, the 



supcrlatticc configuration and the electric field has not 
been addressed earlier. 

To explore it, we present in Fig.|4]thc results of a tradi- 
tional Hubbard chain where all sites are interacting (left 
column) together with the results of a superlattice ge- 
ometry where four non-interacting atoms are placed be- 
tween two interacting atoms (right column). The results 
are shown for three different values of the voltage bias, 
taking a linear bias drop along the ID chain. For the 
chain where all sites are interacting a single energy gap 
only appears at the band centre, while in the superlat- 
tice geometry, depending on the unit cell configuration, 
multiple energy gaps are generated which are clearly vis- 
ible from the ADOS spectra. Therefore, in a superlattice 
geometry, in presence of external electric field associated 
with bias voltage V between two electrodes, zero trans- 
mission (T — 0) energy regions exist, and are separated 
by regions of extended states compared to the traditional 
Hubbard chain, and, it leads to the possibility of getting 
an Mi-like transition at multiple energies. 

The total number of energy sub-bands in a superlat- 
tice geometry for a particular energy range generated in 




FIG. 8: (Color online). Transmission probability T (red color) 
and ADOS (green color) as a function of energy E for a ID 
chain (N = 300) with no disorder (W = 0). The model 
parameters are U = 2, M = 4 and V = 0.2, with the bias 
potential profile taken as the green curve given of Fig. [2] 

the ADOS profile strongly depends on structural details 
i.e., the number M of non- interacting atoms between two 
interacting lattice sites. This is shown in Figs. [S] and [S] 
which show the ADOS and the transmission probabil- 
ity for two models that are identical in all details (see 
caption to Fig. [5]) except that M — 5 in Fig. [5] and 
M = 6 in Fig. [5] These structures show more mobility 
edge phenomena, that is crossovers between fully opaque 
and a transmitting zone, than in the corresponding case 
of Fig. 0|f), suggesting a design concept based on such 
superlattice structures as a switching devices at multiple 
energies. 

The robustness of the observed behavior can be exam- 
ined by its sensitivity to the presence of disorder. Fig- 
ure [7] displays the ADOS spectrum and the total trans- 
mission probability for a ID chain in presence of diagonal 
disorder affected by choosing from a uniform distri- 
bution of width W = 0.5 (-0.25 to +0.5). An aver- 
age over 50 disorder configurations is presented. The re- 
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suiting ADOS and transmission show similar qualitative 
features, with sharp transitions between localized and 
extended spectral regions as seen above for the ordered 
cases. Note that the presence of disorder alone can cause 
state localization. For strong enough disorder almost all 
energy levels get localized, even for such a finite size ID 
chain. In this limit such crossover behavior will no longer 
be observed. 

In the calculations presented so far we have assumed a 
linear drop of the electrostatic potential along the chain. 
Figures |S] and show results obtained for an identical 
chain length with other potential profiles that are char- 
acteristics of stronger screening. We see that the local- 
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FIG. 9: (Color online). Same as Fig. [8] with the electrostatic 
potential profile given by the blue curve in Fig. [2] 

ized region gradually decreases with increasing flatness 
of the potential profile in the interior of the conduct- 
ing bridge. If the potential drop takes place only at the 
chain-to-electrode interfaces i.e., when the potential pro- 
file becomes almost flat along the chain the width of the 
localized region almost vanishes and the metal-insulator 
transition is not observed, as was the case for the zero 
bias limit. 

Finally, we point out that by locating the Fermi energy 
in appropriate places of the sub-bands, the system can 
be used as a p-type or an n-type semiconductor. For ex- 
ample, let us imagine, at absolute zero temperature, the 
Fermi level is fixed in the localized region which is very 
close to the fully transmitting zone (right hand side). In 
this case, the left sub-bands up to the Fermi level are 
completely filled with electrons. Now, if the energy gap 



between the Fermi level, pinned in the localized region, 
and the bottom of the transmitting region (right hand 
side) is small enough for electrons to hop, then the sys- 
tem will behave as an n-type semiconductor. On the 
other hand, by reverting the situation we can generate a 
p-typc semiconductor where electrons hop from a filled 
transmitting zone (valence band) to unoccupied localized 
zone (conduction band) generating holes in the valence 
band. 

IV. CONCLUSION 

To summarize, we have investigated in detail the 
two-terminal finite bias electron transport in a ID su- 
perlattice structure composed of interacting and non- 
interacting atoms. The electron-electron interaction is 
considered in the Hubbard form, and the Hamiltonian is 
solved within a generalized HF scheme. We numerically 
calculate two-terminal transport by using a Green's func- 
tion formalism and analyze the results for some specific 
chain structures characterized by different arrangements 
of the atomic sites in the chain. Our analysis may be 
utilized in designing a tailor made switching device for 
multiple values of Fermi energy (or, more practically, for 
different values of a gating potential). The sensitivity 
of this switching action i.e., metal-to- insulator transition 
and vice versa on the electric field variation has also been 
discussed. Though the results presented in this article 
are worked out at absolute zero temperature limit, the 
results should remain valid even at finite temperatures 
(~ 300 K) since the broadening of the energy levels of 
the superlattice structure due to its coupling with the 
metal electrodes is much higher than that of the thermal 
broadening^. 
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